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SUMMARY 


A three-dimensional guidance scheme for descending from lunar orbit to a 
hovering position was developed and analyzed. The scheme is based on the 
linear theory of perturbations about a nominal reference trajectory and uses 
thrust acceleration and thrust orientation angles as variables in the control 
equations. It allows the preselection of a lunar landing site from a wide 
range of initial conditions and permits guidance within the rigid constraints 
proposed for LEM. The specific control equations are fully developed in the 
report. Guidance capability, the effect on guidance capability of reducing 
the target size, and fuel consumption are considered in detail. An analog 
computer was used for the investigation and the basic results were checked by 
means of a digital computer program. 


INTRODUCTION 


The powered descent from lunar orbit to a hovering position near the 
lunar surface is critical for various reasons. The vehicle will be required 
to reach a preselected hovering point, and restraints on position and velocity 
will be very stringent. The amount of propellant required will influence the 
earth launch requirements and the payload delivered to the lunar surface. 
Excessive fuel consumption during the powered descent would reduce hovering 
time, thus limiting the time allowed for selecting a desirable landing site. 
These are some of the major problems of a lunar descent. 

Various guidance schemes have been investigated for the powered lunar 
descent (e.g., refs. 1-6). A two-dimensional guidance scheme for lunar 
descent was developed in reference 1. It was based on linear perturbation 
theory and two types of reference trajectories, the gravity-turn and constant 
pitch-rate maneuvers. The purpose of this report is to extend the two- 
dimensional gravity- turn study of reference 1 to a full three-dimensional 
analysis. In addition, a different variable, thrust acceleration rather than 
thrust, will be used in the guidance equations because it is more readily 
measured. The control equations thus differ from those of reference 1 and new 
two-dimensional guidance results are presented. The specific equations nec- 
essary for three-dimensional guidance are developed in the appendix. Detailed 
fuel consumption data are presented to show the efficiency of the guidance 
system . 

The effect on guidance capability of reducing the target area and of 
initial errors in vehicle flight-path angle is also shown. 



The performance of the guidance system was investigated by simulating the 
system and mechanizing the equations of motion on an analog computer. The 
basic results were checked using a digital computer program. 

The authors are indebted to Dr. Richard Rosenbaum of Lockheed Missiles 
an d Space Company for pointing out the advantages of using thrust acceleration 
as a control variable . 


NOTATION 


Ap thrust acceleration, g 

FiJ linear theory gain for the x state variable used to determine the 

magnitude of the control variable, u, dimensions of u/x 

g e earth surface gravity 

h altitude above lunar surface , m 

I S p specific impulse, sec 

m mass, kg 

r r m + h, km 

r m radius of moon, km 

T thrust , N 

V total velocity, m/sec 

X r m ¥, ran § e -’ 

XqiQ. range to go , km 

Y r m A, lateral range , km 

7 flight -path angle (sketch (b)), deg 

S( ) difference between actual and reference value, ( ) - ( ) r 

£ heading angle ( sketch (b ) ) , rad 

9 thrust angle (sketch (c)), rad. 

A lateral range angle (sketch (a)), rad 

^x GrO variable 
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\i product of universal gravitational constant and mass of planet , m 3 /sec 2 

¥ range angle (sketch (a)), rad 

^ thrust angle (sketch (c)), rad 

(*) derivative with respect to time 

( )' derivative with respect to velocity 

Subscripts 

a actual 

f final 

i initial 

j *k summing indices 

r reference 


LUNAR LANDING APPROACH 


Figure 1 is a schematic diagram of a lunar landing maneuver. It is 
assumed that a manned space vehicle enters a l6o km (100 statute miles) circu- 
lar orbit about the moon from a direct earth-moon trajectory. The vehicle 
then performs a transfer maneuver, which places it in an elliptical lunar 
orbit with a perilune of 22,900 m (75*000 ft) and makes a powered descent 
from perilune. 

The Hohmann transfer maneuver is assumed here; however, any transfer 
maneuver that will place the vehicle in a lunar orbit with a perilune of 
22,900 m (75*000 ft) could be used. The guidance scheme in this report is 
developed for the descent from perilune and does not include the transfer 
maneuver . 


REFERENCE TRAJECTORY 


The nominal or reference trajectory chosen as a basis for developing the 
guidance system is the same constant-thrust gravity-turn descent trajectory 
that was considered in reference 1. It begins at perilune of the elliptic 
transfer orbit and terminates at the hover altitude of 300 m (l,000 ft) . The 
initial altitude is 22,900 m (75*000 ft) and the trajectory traverses approx- 
imately 300 km (190 statute miles) of lunar surface. Pertinent quantities of 
the reference trajectory are shown in figure 2. 
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The gravity-turn reference trajectory was chosen primarily because it is 
convenient and requires less than 5 percent more characteristic velocity than 
an optimum trajectory beginning at the same perilune altitude (ref. 7) • 


PERFORMANCE CRITERIA 


In this investigation the target area was assumed to be 300 kni (190 miles) 
from the reference trajectory perilune. The dimensions of the target area are 
the same as those used in reference 1, that is, altitude, ±150 m (±500 ft), and 
range, ±1,500 m (±5,000 ft). An additional dimension of ±1,500 m (±5,000 ft) 
is assumed for crossrange. The center of the target area is at the hover alti- 
tude of 300 m (l,000 ft) . Performance capabilities, presented elsewhere in 
this report, are based on the ability of the vehicle to reach this target area 
with less than 3 mps (10 fps) vertical velocity and less than 3 nips (10 fps) 
total velocity. Performance capabilities based on target areas that are l/2 
and l/4 of the aforementioned target area are also presented and compared with 
the capabilities based on the nominal target area. The comparison gives some 
indication of the effect of target area dimensions on the guidance system 
capability. 

No performance criteria were selected for fuel consumption; fuel consump- 
tion data for the trajectories investigated are presented simply to indicate 
the fuel necessary to attain the guidance capabilities. 

INITIAL CONDITIONS 


In this investigation initial perilune altitudes were considered up to 
±9,200 m (±35,000 ft) from the reference and initial ranges, 240 to 900 km 
(150 to 560 miles) from the target. The reference trajectory range is approx- 
imately 300 km (190 statute miles)-. Initial lateral displacements up to 
l60 km (100 statute miles) outside the plane of the reference trajectory were 
also considered. The expected errors in altitude, range, and crossrange for a 
LEM mission are considerably smaller than the spread of initial conditions 
considered here. 


GUIDANCE EQUATION 


The concept of guidance about a nominal or reference trajectory is based 
on perturbation theory, that is, the analysis of conditions in a limited 
neighborhood of a nominal trajectory. The basic linear perturbation theory 
guidance equation is developed in reference 1. To control altitude above the 
lunar surface, h, range, X, and lateral range, Y, at the hover point, the 
quantities, thrust acceleration, Aj,, thrust pitch angle, 0, and thrust yaw 
angle, i|r, are varied. Thrust acceleration rather than thrust is the control 
variable because it is more easily measured, the number of terms in the guid- 
ance equations are reduced, and deviations in thrust and mass from their 
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nominal values are accounted for by one readily available measurement. Thus, 
fewer guidance gains are required and less onboard computer storage is 
required. 

The specific guidance equations used in this investigation are developed 
in the appendix. They are: 


Aj(v) = A* (V) + f£ T (V)S1i + F^ t (V)Sx + Fy T (V)57^ 
e(v) = e r (v) + F®(v)sh + f®(v)b x + f®(v)&7 > 

t(v) = i|r r (V) + Fy(V)BY + fJ?(V)B£ 


( 1 ) 


where total velocity, V, is the independent variable, 5( ) is the difference 
between the actual and reference value, ( ) a - ( ) r > and F are the guidance 
gains defined in the appendix. 


Total velocity, V, rather than time, was chosen as the independent 
variable to reduce the number of terms in equations (l) by one, thus simplify- 
ing the equations and reducing information storage requirements. 


A flow diagram of the guidance system operation is given in figure 3 in 
terms of a general guidance equation. Note that the reference state variables, 
reference control variables, and guidance gains are stored as functions of the 
independent variable velocity. Control increments are computed on the basis 
of differences between the actual trajectory and the reference trajectory at a 
given velocity. 

The guidance equations with thrust acceleration, Ap, used as a control 
variable allow the effect of perturbations in thrust and mass from their 
nominal values to be combined into a single term AA p. 


GUIDANCE GAINS 


The calculated linear theory guidance gains associated with equations (l) 
are shown in figure 4 . The performance capability was investigated by simu- 
lating the guidance system on an analog computer. Because of the limited 
storage capacity of this computer, the guidance gains were stored for 17 
velocity points with linear interpolation between points. To insure accu- 
racy in storing the gains within the scaling and storage restrictions of the 
analog equipment, the gains were limited as indicated by the dashed lines of 
figure 4. 

An important result of reference 1 was that the guidance capability could 
be greatly increased by weighting the guidance gains empirically. In the 
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present study the gains in figure 4 were increased and decreased individually, 
and it was found that increasing the gain F^ by a factor of 2.5 and F A by 
a factor of 1.5 substantially improved the guidance capability. It should be 
emphasized that these factors are not optimum in any way. They are merely two 
of the more effective values of the guidance gains discovered by increasing 
and decreasing the gains individually. The adjusted gains were limited at the 
same value of velocity as were the gains in figure 4, and no attempt was made 
to adjust more than one gain at a time. There are various methods of adjust- 
ing the gains; one simple method would be to change the value of velocity at 
which the gains were limited. Conceivably there are other combinations of 
adjustments which would similarly increase capability but no attempt was made 
to find any optimum set of guidance gains . The intent is simply to show that 
the guidance capability can be improved in a manner similar to that 
demonstrated in reference 1 . 

An IBM 7094 digital computer program was used to check the analog computer 
results. Many more than 17 points were used, making it possible to ascertain 
the validity of the simpler storage of the analog computer . Since the digital 
n-nrl analog results compared quite well, the analog storage is considered valid. 

An attempt was made to further simplify the guidance gains by approximat- 
ing each gain shown in figure 4 with a single straight line . Several differ- 
ent combinations of single straight line approximations were investigated but 
in each case the guidance capability was greatly reduced in comparison with 
the linear theory capability. This gives some indication of the importance of 
the shape of the guidance gain curves. 


VEHICLE CHARACTERISTICS AM) EQUATIONS OF MOTION 


If the vehicle is assumed to be a point mass, and the moon to be 
spherical and nonrotating, the resulting equations of motion are: 

h = V sin 7 

7 = 4 1 sin 9 cos + V cos 7 _ i JL CO s 7 
mV r V r^ 


V = ?. cos d cos j _ iL sin 7 


m 


i = 1 s - -4 n . -ty - ■ - — tan A cos 7 cos £ 

mV cos 7 r 


Y = r A = — V cos 7 sin £ 


m 


X = r m T 


• _ £m V cos 7 cos C 


cos A 
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The mass rate is given "by 
m = mi - J*(T/g e I sp )at, where I sp 
equals 425 seconds . The geometry for 
the equations of motion is shown in 
sketches (a), (h) , and (c). Sketch (a) 
presents the axis system for the tra- 
jectory equations. The coordinate sys- 
tem is centered at the center of the 
moon. Sketch (b) defines flight-path 
angle 7 and heading angle £ . Thrust 
angle orientation is shown in sketch (c). 
Thrust pitch angle Q is the angle 
between the velocity vector and the 
projection of the thrust vector in the 
vertical plane . The thrust yaw angle 
t is the angle between the thrust 
vector and the vertical plane . 

The vehicle is assumed to have a 
gimballed, variable -thrust rocket 
engine, ho throttle ratio was imposed 
but the minimum thrust was limited to 
1.11x10 4 N (2,500 lb) which was one- 
tenth the reference thrust. The guid- 
ance system actually required a maximum 
thrust of 156,000 N (35,000 lb) to 
achieve the full three-dimensional capa- 
bility shown elsewhere in this report. 

A maximum thrust of 142,000 N (32,000 
lb) was required for the two-dimensional 
capability. Two-dimensional capability 
for ranges longer than the reference 
range required maximum thrust to be less 
than 111,000 N ( 25,000 lb). 


RESULTS 


Two-Dimensional Guidance Capability 

The two-dimensional guidance 
capability of the guidance system when 
the descent maneuver is initiated at 
perilune of an elliptic Hohmann trajec- 
tory is illustrated in figure 5(a) . 

The capability is presented in terms of 
the initial altitude and range limits 
from which it is possible to meet the 
altitude and range constraints of the 
aforementioned target area. Initial 
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values of altitude and velocity correspond to the Hohmann trajectory with 
perilune at that point. Since the descent was always initiated at perilune, 
the vehicle's initial flight-path angle and vertical velocity are zero. 

The smaller capability indicated by the shaded area was obtained from the 
linear theory guidance gains of figure 4. The increased capability indicated 
by the boundaries labeled "limits of investigation," which include the shaded 
area, was obtained by increasing the gain Fj^ by a factor of 2.5* 

The heavy black lines of figure 5(a) indicate limits of capability 
defined by the investigation. The target constraints which were exceeded are 
also indicated. 

For the guidance capability in figure 5(a) numerous initial points were 
considered at three different perilune altitudes; the reference value and the 
two altitudes corresponding to the altitude limits of investigation, namely, 
12,200 m (40,000 ft) and 33,500 m (110,000 ft). Initial range-to-go values 
were varied at intervals of approximately 40 km (25 statute miles) across the 
region of capability. A brief check was made of the results for trajectories 
whose perilune altitude was intermediate to the three discussed above and also 
of trajectories whose initial range to go was intermediate to the intervals of 
40 km. Results for these intermediate trajectories correlated well with the 
basic data. The consistency of results for this large amount of data indicates 
that guidance to the target area could be achieved from any point within the 
region of capability. 

Figure 5(b) is a comparison of the guidance capability using the linear 
theory gains with capability as presented in reference 1 using linear theory 
gains. As Indicated, the capability is generally greater for the guidance 
scheme presented in this report. The improved capability is a result of 
taking into account vehicle mass perturbations in the development of the 
guidance equations. The guidance results from using empirically adjusted 
gains are not compared directly as it would be extremely difficult to determine 
which set of guidance equations would yield greater extended capability. The 
intent here is to verify that the capability can be extended in the manner 
shown in reference 1 . 

Figure 5(c) compares the guidance capability for the nominal target with 
that for two targets, the dimensions of which are one -half and one -fourth of 
the nominal dimensions. As indicated, the capability is substantially reduced 
when guiding to a smaller target. The data in figure 5(c) are for the linear 
theory gains . 

In figure 5(&) the same comparison is made for capability using the 
empirically adjusted gains. In this case the capability is not altered by 
target size . 

The two-dimensional guidance capability of the system when the descent is 
Initiated at a point other than perilune is presented in figure 6 in terms of 
the initial altitude and range limits from which it is possible to meet the 
altitude and range constraints of the target area. Since the descent is not 
initiated at perilune, initial flight -path angle and vertical velocity are 
not zero. For the purposes of this phase of the investigation it was assumed 
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that the Hohmann transfer is always initiated at the same apolune point l60 km 
(100 statute miles) above the lunar surface. Therefore, initial conditions are 
determined by the Hohmann trajectory passing through that point. Values of 
flight-path angle along three of the Hohmann trajectories considered are given 
in figure 7 . Initial flight-path angles ranging from -0.8° to 0,1° were con- 
sidered. The corresponding vertical velocities ranged from -23 m/sec to 
6 m/sec (-75 ft/sec to 20 ft/sec). 

The smaller capability indicated by the shaded area resulted from using 
the calculated gains of figure 4 and the increased capability indicated by the 
boundaries labeled “limits of the investigation/ 1 which include the shaded 
area, was again obtained by increasing the gain F^ by a factor of 2.5. Again 
the heavy black lines indicate defined limits of capability and the constraints 
which define the limits are also indicated. 

The guidance capability presented in figure 6 was determined in a manner 
similar to that used for figure 5(a). Descents were initiated at numerous 
points along three trajectories whose perilunes were located at the reference 
range from the target. The perilune altitudes were 12,200 m (40,000 ft), 

22,900 m (75,000 ft), and 30,500 m (100,000 ft). The interval between values 
of initial range to go was again approximately 40 km (25 statute miles). Again 
the consistency of the data indicates the vehicle can reach the target area 
from any point within the region of capability. 

In both figure 5(a) and figure 6, the only limit of capability when the 
gain 2.5 F^ is used is for shorter ranged trajectories than the reference. 

This limit was defined by failure to meet range constraints. Elsewhere the 
capability extends to the limits of the investigation . 

Three-dimensional guidance capability .- The results discussed thus far 
are for two-dimensional trajectories. How consider trajectories which have an 
initial lateral displacement out of the plane of the reference trajectory. 
Descent is initiated at perilune. The initial conditions of figures 8(a), 

8(b), and 8(c) are identical to those of figure 5(a) with an added initial 
lateral displacement. Figures 8(a), 8(b), and 8(c) present the lateral guid- 
ance capability for three different perilune altitudes in terms of the initial 
lateral displacement and range limits from which it is possible to guide to 
the specified target area. The smaller capability, indicated by the dashed 

lines, was obtained using the values of F^ and F^ given in figure 4. The 
wider capability for each altitude, indicated by the solid lines, was obtained 
by increasing F^ by a factor of 1 . 5 . This factor was one of the more effec- 

* \lr Mf 

tive values discovered by increasing and decreasing Fy and F£ separately. 

Other combinations of gain adjustments might give equally good increased 
capability but, again, no attempt was made to obtain an optimum set of gains. 

In all three-dimensional work, the gain 2-5 was also used. 

The constraints which were exceeded in defining the limits of capability- 
are indicated in figure 8. 
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A comparison of figures 8(a), 8(b), and 8(c) shows that lateral capability 
is greater at lower altitudes for short ranged trajectories and greater at 
higher altitudes for longer ranged trajectories when the gain 1.5 F^ is used. 
When the linear values of is used, the capability is generally larger for 

higher altitudes. The maximum lateral capability is 160 km (100 statute miles) 
initial displacement. This is the lateral capability for a descent whose 
perilune is at 33, 500 m (ll0,000 ft) altitude located 900 km (560 statute miles) 
from the target. In general, the system has a large lateral capability with 
linear theory gains which may be increased with adjusted gains. It will be 
pointed out in the next section that fuel consumption is greatly Increased for 
trajectories having initial lateral displacements. This increased fuel 
consumption may limit the useful lateral capability. 

An important conclusion of the study of lateral range guidance is that 
the addition of the third dimension to be controlled does not reduce in any 
way the two-dimensional guidance capability. 

FUEL REQUIREMENTS 


Figures 9(a), 9(b), and 9(c) present the additional characteristic 
velocity, AV C , required for the lateral range capabilities with the adjusted 
gain 1.5 f| in figure 8. Additional velocity, AV C , is defined as V c of the 
actual trajectory minus V c of the reference trajectory, AV C = V Ca - V Cr , 
where V c = g e I S p 2n(mp/mp) . The curves shown represent AV C as a function 
of initial lateral range, Yj_, and initial range to go, XpQ^. The figures indi- 
cate that the guidance system requires little AV C if the vehicle is not 
required to change orbital plane. The maximum A V c required for the two- 
dimensional trajectories of figure 5 is only 21 m (70 ft) per second. The sys- 
tem can cope with initial lateral displacements of 24 to 40 km (15 to 25 
statute miles) with only a small additional AV C . However, as initial lateral 
displacement increases, the characteristic velocity, which is a measure of 
fuel consumption, becomes an important consideration and may limit the lateral 
capability. The large increases in AV C with initial lateral displacements 
result from the necessity to change the plane of the trajectory. 


CONCLUDING REMARKS 

A three-dimensional guidance scheme for the descent from lunar orbit to 
hover has been investigated. The guidance scheme was developed by applying 
linear perturbation theory to a single nominal or reference trajectory. It is 
shown that the scheme has guidance capability beyond that required for lunar 
descent of the type proposed for the LEM. Thrust acceleration, At, a readily 
measured quantity, and thrust orientation angles are used as control variables. 
With this scheme a hover point near the lunar surface can be preselected from 
a wide range of initial conditions; empirical weighting of the computed guid- 
ance coefficients would compensate for initial ranges of 240 to 900 km (150 to 

560 statute miles) from the target, from altitudes of 12,200 to 33*500 m 

(40,000 to 110,000 ft), and initial lateral displacements up to l60 km 

(100 statute miles) outside the plane of the reference trajectory. 
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Range and altitude guidance is accomplished with only a nominal increase in 
fuel consumption hut lateral displacements require more . 


An important result of this study is that the addition of a third 
dimension to the guidance does not reduce the two-dimensional capability. 

Empirically adjusted gains permit guidance to a much smaller target area 
with only a small decrease in guidance capability. 

Since only a single nominal trajectory is used, little onboard computer 
storage capacity is required. 

Ames Research Center 

National Aeronautics and Space Administration 
Moffett Field, Calif., Sept. 7> 19^5 


11 



APPENDIX 


DEVELOPMENT OF SPECIF IG CONTROL EQUATION 


It is desired to develop a terminal guidance system based 
tion of linear perturbation theory to the reference trajectory, 
perturbation equations which describe the motion of a body near 
jectory are the total differentials of the equations of motion, 
the form: 


on the applica- 
The linear 
a nominal tra- 
They are of 


dp = Sh + ^ 5x + ~ &7 + M SV + ^ Sy + ^ 
oh 9x 07 oV ay 


+ 6A T + M 50 + M St 
SAqi S0 <3^ 


(Al) 


or 


5p = A-j^Sh + A i2 &X + A^3&7 + Aj_ 4 5V + A^ 5 8y + A i6 &£ + bp-,SAp + b^ 2 50 + b i3 S\|r 

(A2) 

where 

P = A, x, 7, V, Y, i 

and 

i = 1^ 2^ 3; 5; ^ 


The control variables are thrust acceleration, Arp, thrust pitch angle , 

By and thrust yaw angle , The use of thrust acceleration as a control vari- 
able simplifies the control equations by allowing the effect of perturbations 
in mass and thrust to be combined into a single term, AA«p. 


Total velocity, V, is the independent variable in the control equations. 
The equations of motion are written thus: 


or 


P = 


dp = d£ dV 
dt dV dt 


= p T V 


P 1 


= P 
V 


(A3) 


where p = h, X, 7 , A, £ and the prime refers to differentiation with respect 
to total velocity., V. 


The linear perturbation equations now become: 


, , dp 1 

dp ’ 5T 


Sh + 


dp 1 

dx 


sx + 

^7 


87 + 


8 p' 

dF 


6 Y + 


s? + 5At + 5£ > + hi 

d£ dArp de dt 


6 f 

(A4) 
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or 


Bp* = A kl Sh + A k2 8X + A k3 &7 + A k4 SY + + B^SAj + + Bj^St 

(A5) 

where 

p = h, X, 7> Y, £ 

and 

k = 1; 2, 3, b } 5 


The partial derivatives are evaluated for the constant -thrust gravity- 
turn reference trajectory. By definition, 9 = 180° on the reference trajectory 
and, since a planar trajectory is assumed, A = £ = = 0 . The nonzero coeffi- 

cients of equation (A5) are found to he: 


li 


*■13 


8h' -li 2p . ^ 

- — = t-v — v sin 7 

8h V 2 -3 


8h' 

87 


r 

V cos 7 
V 


M- 


h 


+ — =r cos 7 


l 21 


^23 


'SI 


_ SX' _ -r m V cos 7 


8h 

8X' 

87 

87' 

8h~ 


r 2 tf 

r m V sin 7 


rV 

-V cos 7 


% sin 7 \ 
r s ^2 

+ -4 cos 7 4 


V 


2(i 

+ v?3 cos y 


. 2p 

7 -5- cos 7 
r 2 


V 


a = £!2Ll 
Ass 87 


sin 7 + ^ ^ sin 7 ^5 cos 77 

1 — T — + — 


M5 


= 8Yj_ = fm V cos 7 


8£ 


V 


A, 


54 


Bn = 


B 2 i - 


-V cos 7 


8Y 

8h' 

8Aj v 


rV 


.2 


8x» 

8Ap 


X 


T>„ - ^Zl - _L 

B31 ‘ sat ’ v 2 


13 



T 


b 32 = |21L = 

8e mW 


B 


53 


<&: 

ch|r mV cos TV 


In matrix form the linear perturbation equations are: 


8h* 


An 

0 

Al3 

0 

0 


5h 


Bn 

0 

0 


6X' 


A 2 i 

0 

A 23 

0 

0 


5X 


B 21 

0 

0 


5Ap 

57 ' 

= 

A 31 

0 

A 33 

0 

0 


67 

+ 

B 31 

B 32 

0 


50 

SY' 


0 

0 

0 

0 

in 

<J 


SY 


0 

0 

0 


_s t 

l51 


0 

0 

0 

A 54 

°J 




0 

0 

B 53 



(a6) 


They are linear differential equations with variable coefficients. It 
should be noted that the last two equations are uncoupled from the others, 
indicating that, to first order, the heading angle £ and lateral displacement 
Y are independent of range, altitude, and flight -path angle. This is only 
true for a planar reference trajectory. 


Values of 6 h, SX, S 7 , 6 Y , and S£ at the final velocity, V^ , can be 
determined, if their values are known at some initial velocity, Vq, by setting 
SAp = 50 = SV = 0 . The effect of a control variable on the final values of the 
state variables can be found by solving the perturbation equations with all 
other variables initially equal to zero. 


This method of computing control increments must be repeated for every 
new Vq because of the variable coefficient equations. It becomes extremely 
cumbersome to calculate control increments at every point along the trajectory. 
This may be avoided if the method of adjoint functions is used (see ref. 8). 

One solution of the equations adjoint to the linear perturbation equations 
yields all the required information. 


In matrix form the equations adjoint to the linear perturbation equa- 
tions (A6) are: 


K 


An 

A21 

Asi 

0 

0 

*h 



0 

0 

0 

0 

0 

*x 


= - 

A 13 

A23 

A33 

0 

0 




0 

0 

0 

0 

A54 

A Y 

_ a L 


L° 

0 

0 

A45 

°J 

aJ 


(AT) 
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Equations (A 6 ) may be -written 


13 3 

x i - X *kj x J * I V) p j 


(A 8 ) 


J=i 


0 =1 


and equations (A7) may be written 


Nl + X a j 


j =1 


... A, = 0 

jk J 


(A9) 


Now multiply (A 8 ) by A^. and (A9) by X x ; then add and sum over k: 

5 _ 5 ^ 5 5 3 

X ^ A k X k + X k A k) + X Yj ^ X k a jk A j " x j a kj A k) = X X Akbk J p j ^ A1 °^ 


k=i 


j=i k=i 


k=i j=i 


The second term in (A10) is zero since j and k are summing indices over the 
same range. Therefore (A10) becomes 


5 3 


X ( A k x k + x k\) = X X Akb 


_ _ kj p j 

k=i k=i j=i 

The left-hand side is a perfect differential so write (All) 

5 5 3 

§ X ( w - X X Vkjpj 

k=i k=i j=i 

Integrating (All) from V = V 1 to V = Yf gives 

5 3 


(All) 


X \(V f ) x k(V f ) = ^ x k( v r) x k( v i) + f f \ X A k( v ) b kj( V )Pj( V ) dV ( A12 ) 

i i, — i V_ i, — , ■ 


k=i 


k=i 


i k=i 3=1 


Equation (A12 ) , Bliss' fundamental formula, is the basic equation for control 
about a reference trajectory (see ref. 9 ) • 

Substituting values from (A7) and (A 8 ) into (A 12 ) gives the following 
equation: 

[AftSh + A X 5X + Ay 5y + A y 5Y + A^5^] y=Vf 

= [A^Sh + A^SX + Ay 67 + A Y 5Y + A£S£] y _ y 

r v f 

+ J (b 32 Ay50 + b 41 A Y 8 A T + b 63 A^5t)dV (A 13 ) 
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To control the final altitude, h, range, X, and lateral displacement, Y, 
we must solve equation (AI3) for the final errors 5hp, SXf, and SYp at 
V = Yf . Note that the left-hand side of (AI3) equals 5hf (if Ah = 1 and all 
other A' s are zero). Solving (AI3) in this manner for Shp , 5Xf, and 5Yf 
yields the following expressions: 


Shf = 


SX-F = 


8Yf = 


a|V + A^Sy 


+ 


iy=Vi 0 V 


rVf h pVf 2i 

J (Bi 1 A li + B3iA 7 )8A T dV+ J B 32 A 7 50 dV 


A^8h + 5X +A^67 




f v Y P X V 

+ / (Bh\ + B21 + B 3iA^, )dV + / B32A 80 dV > 

J V =Vi J Vx n J Vi / 


A^SY + A^5(; 




J V=Vi 


B 53 A£&\|r dV 


J 


(Al^) 


where the superscripts on the A' s refer to the final error being solved for. 
These equations give the terminal values of the variables to be controlled as 
functions of initial deviations and of the control variables. 


Given a desired final value of the state variables to be controlled, 
there are any number of control variable functions which can be used to attain 
the desired final values. In particular, there are constant values of the 
control variables over the interval V < Vi < Vp which will result in the 
desired final values of the controlled state variables. The linearity of the 
equations allows the addition of the separate conditions to determine the con- 
stant values of the control variables that will eliminate off -design conditions 
in the controlled state variables at the final velocity. The guidance system 
then calls for control increments which, if held constant, will just eliminate 
errors in the controlled state variables at the final velocity. 

Equations (Al4) will be solved for the constant control variables AAj 1 , 

A 0, and A\|r. A shorthand notation will be introduced for the sums and integrals 
of equations (Al^) as follows: 


E SH = A^h + ^7 S A 


Egx = A^8h + SX + A^6 7 


Jy=Vi 

X, 

Y 


ly=Vi 


E sx - 

X TH = 


A|5Y 


A^5£ 


J V=Vi 



+ B 



dV 


l6 



Equations (Al4-) may now be written: 

8 hf = E SH + ^T^TH + S0I 0H 
5X-f = ^SX + ^Ajplrj-r^- + &Q3-QX 


8Y f = Egy + 


(A15) 


If Sj^ and 5>Yf are to be zero , equation (A15) must be solved for the control 
variables: 

= ( E sx - SXf) 1 !^ ~ e SH I 0X 
IthIqX - ItxIqH 


80 = 


8 i|r = 


e SH I TX “ ( E SX " SXf) 1 ^ 


I TH I 0X - I TX I 0H 


~ E SY 

I 'l r Y 


(Al6) 


When equations (Al6) are expanded and like terms are collected, equations (Al6) 
may be simplified as follows: 


SAjji = 


50 = ^ 


A h I 0H " A h I 0X 8h " I 0H sx TG + A 7 X 0H " ^ax 8 ? 




X h T 

A h ± TX 




ITHI0X - I TX I 0H 
5h + IijijjbXrjTQ. + 


A 7-'-TX “ A 7 E TH | 8 7 


lTH I 0X " I TX I 0H 


&t = - 


r 7^SY+A Y 


I tY 




(A 17 ) 
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Illllllllllllll 111 


Equations (A17) will "be used to develop the guidance system. In equations (A. 17) 
let 


F. 


t A T ^OE " 


F 


h 

Arp 

r x 

7 


D 


-SflB 

D 


*\ X-r- -\ll-T- 

A 7 J ’9H ~ A 7 x 9X 
D 

_ ^h^x " 

f h = 5 ^ 


v e - 
i x - 


^TH 

D 


■nd _ ^i^x - ^th 
7 D 

f* = 'h. 

i i tY 

-vY 

F t = 

^ I \(ry 


where 


D = Ipn-Lex " I TX I 0H 

Substituting in equations (A17) gives: 

Atji Aip Arp 

XAji — F^ &li + F-^ ^7 

A 6 = F®6h + F^SX^ + F®&7 


\[r \Lr , 

A = FySY + FpSf; 


(A18) 


These equations are used to compute control variable corrections and; 
when added to reference trajectory values of the control variables; give the 
following control equations: 
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(A19) 


Arp A rp Aj 1 

ApCv) = A^V) +F h (V)5h +F X (V)SX TG + F y (V)5/ 1 

e(V) = 0 r (v) +F®(¥)5h +f|(y)5X tg +F®(v)&7 > 

*(V) = t r (v) + f|(v)BY + F^(v)&£ j 

These are the specific control equations used in this investigation. 
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Figure 1.- Linar landing maneuver. 
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Figure 2.- Continued. 




(c) Characteristic velocity and thrust acceleration. 


Figure 2.- Concluded. 




Figure 3.- Terminal control system; where U 


Au p (V)=2F$JP(V)8Xm(V) 
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Xm r (V) 

I 

8Xm( V)=Xm a (V)-Xm r (V) 
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Figure ^ . - Linear theory guidance gains. 
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(a) Descent initiated at perilune . 
Figure 5-- Guidance capability. 
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(b) Comparison of guidance capability based on linear theory gains. 


Figure 5 • ~ Continued . 
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(d) Comparison of guidance capability using empirically adjusted gains for three different targets. 


Figure 5 • - Concluded. 
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Figure 6.- Guidance capability - descent initiated at point other than perilune . 
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Figure 7- Initial flight -path, angles investigated for data of figure 6 
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(a) Perilune at 33,500 m (110,000 ft). 
Figure 8.- Lateral guidance capability. 




(b) Perilune at 22,900 m (75^000 ft). 
Figure 8.- Continued. 
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(c) Perilune at 12,200 m (40,000 ft). 
Figure 8.- Concluded. 
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